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Abstract. The Earth matter effects on supernova (SN) neutrinos can be identified 
at a single detector through peaks in the Fourier transform of their "inverse energy" 
spectrum. The positions of these peaks are independent of the SN models and therefore 
the peaks can be used as a robust signature of the Earth matter effects, which in turn 
can distinguish between different neutrino mixing scenarios. Whereas only one genuine 
peak is observable when the neutrinos traverse only the Earth mantle, traversing also 
the core gives rise to multiple peaks. We calculate the strengths and positions of these 
peaks analytically and explore their features at a large scintillation detector as well as 
• at a megaton water Cherenkov detector through Monte Carlo simulations. We propose 

a simple algorithm to identify the peaks in the actual data and quantify the chances 
O ' of a peak identification as a function of the location of the SN in the sky. 
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1. Introduction 

The neutrino spectra that arrive at the Earth from a core collapse supernova (SN) 
have information about the neutrino masses and mixings encoded in them. The ~ 20 
neutrinos observed from SN 1987A were extensively used trying to constrain the solar 
neutrino parameters as well as #13 and the neutrino mass hierarchy PJ El El Ej The 
determination of neutrino parameters depends crucially on our understanding of the 
primary neutrino fluxes produced inside the SN. In spite of large uncertainties on these 
primary fluxes, some of their robust features may be exploited to identify the type of 
neutrino mass hierarchy and put bounds on the mixing of v e in the "third" neutrino 
mass eigenstate [SUSHI- 

When neutrinos pass through the Earth before arriving at the detector, their 
spectra may get modified due to the Earth matter effects. The presence or absence 
of these effects can distinguish between different neutrino mixing scenarios jHj. The 
comparison of neutrino spectra at two different detectors can clearly give signatures of 
the matter effects, which can be used not only for the determination of the neutrino 
parameters jUJ CD], but also to extract information about the density structure of 
the Earth core llj. The measurements of the Cherenkov glow at IceCube may also 
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be combined with the signal at a water Cherenkov detector like Super- or Hyper- 
Kamiokande to identify the Earth effects [T2] . 

It is also possible to ascertain the presence of these matter effects using the signal 
at a single detector. It has recently been pointed out [T3| that the Earth matter effects 
on supernova neutrinos traversing the Earth mantle give rise to a specific frequency 
in the "inverse energy" spectrum of these neutrinos. This frequency, which may be 
identified through the Fourier transform of the inverse energy spectrum, is independent 
of the initial neutrino fluxes and spectral shapes. Therefore, its identification serves as 
a model independent signature of the Earth matter effects on SN neutrinos, which in 
turn can distinguish between different scenarios of neutrino masses and mixings. 

If the SN neutrinos reach the detector "from below," they have to travel through 
the Earth matter. If the nadir angle of the SN direction at the detector is less than 
33°, the neutrino path crosses the Earth core. An investigation of the effect of the core 
on the observed neutrino spectra is therefore necessary for a complete understanding 
of the Earth effects. As we shall see later in this paper, the passage through the core 
increases the chances of the identification of the Earth effects. The core density is almost 
twice the mantle density and this sudden density jump gives rise to new features in the 
spectra. 

When the SN neutrinos traverse the Earth core in addition to the mantle, one does 
not get a single specific frequency as in the mantle-only case, indeed as many as seven 
distinct frequencies are present in the inverse energy spectrum. We study the strengths 
of these frequency components analytically and show that three of these frequencies 
dominate. These three frequencies can also be clearly observed in the Fourier transform 
of the inverse energy spectrum when averaged over many simulated SN neutrino signals. 

Although it is difficult to isolate these frequencies individually from the background 
fluctuations from a single SN burst, we suggest a procedure that can identify the presence 
of these frequency components in a sizeable fraction of cases. Certain characteristics 
of the distribution of the frequency components in the background fluctuations are 
identified and used to reject the null hypothesis of the absence of Earth effects. 

We quantify the efficiency of this algorithm by simulating the SN neutrino signal at 
a large scintillation detector like LENA ^3] and at a megaton water Cherenkov detector 
like Hyper-Kamiokande. Whereas the scintillation detector has the advantage of a much 
better energy resolution, this is compensated in part by the larger number of events in 
a megaton water Cherenkov detector. 

This paper is organized as follows. In Sec. 121 we discuss the positions and 
the strengths of the frequencies that characterize the "inverse-energy" spectra of the 
neutrinos crossing the Earth mantle as well as the core. In Sec. |3J we simulate the 
SN neutrino spectra at the detectors and study the features of the peaks with the 
background fluctuations averaged out. In Sec. IH we introduce a method to identify the 
peaks in the presence of background fluctuations and make a quantitative estimation of 
the probability of peak identification as a function of the location of the SN in the sky. 
In Sec. we summarize the results. 
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2. Frequencies contributed by Earth effects 

2.1. Mixing scenarios and Earth effects 

The neutrino detectors, apart from a heavy-water detector like SNO, can give detailed 
spectral information only about the v e flux. We shall therefore concentrate on the v e 
spectrum in this paper. In the presence of flavor oscillations a u e detector actually 
observes the flux 

F°(E) , (1) 

where Ff and Ff stand for the initial and detected flux of z/j respectively, and p D (E) is 
the survival probability of a v e with energy E after propagation through the SN mantle 
and perhaps part of the Earth before reaching the detector. The bulk of the v e are 
observed through the inverse beta decay reaction v e p — > ne + . The cross section o of 
this reaction is proportional to E 2 , making the spectrum of neutrinos observed at the 
detector N(E) oc aF s D oc E 2 FP. 

In the absence of Earth effects, the dependence of the survival probability on E 
is very weak. A significant modification of p D due to the Earth effects takes place 
only when the neutrino mass hierarchy is normal, i.e., mi < m-i < 777,3, or when the u e 
component of the third mass eigenstate is restricted to \U e ^\ 2 = sin 2 13 < 1CT 3 . Here 
^3 is the neutrino mass eigenstate that has the smallest v e admixture. The identification 
of the Earth effects can then rule out the "null hypothesis" of an inverted hierarchy and 
I t/ e 3 1 2 > 10~ 3 , thus excluding a large chunk of the neutrino mixing parameter space. In 
the language of Table ^ the Earth effects can be present in scenarios A and C whereas 
they are absent in scenario B. Case B is thus the null hypothesis. 



Case 


Hierarchy 


sin 2 0i3 


Earth effects 


A 


Normal 


Z lo -3 


Yes 


B 


Inverted 




No 


C 


Any 


^ io~ 3 


Yes 



F°(E)=p D (E)F° B (E) 



(E) 



Table 1. The presence of Earth effects in different neutrino mixing scenarios. 



Let us consider those scenarios where the mass hierarchy and the value of G13 are 
such that the Earth effects appear for v e . In all of these cases, v e produced in the SN 
core travel through the interstellar space and arrive at the Earth as V\. The oscillations 
inside the Earth are essentially V\-v-i oscillations 5\ so that we need to solve a 2 x 2 
mixing problem. 



2.2. Passage through only the mantle 



When the antineutrinos pass only through the mantle which has roughly a constant 
density, the survival probability p D is given by 

2 
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R(-e m )<b(<p m )R(6 m - e 



vi. 



11 



(2) 
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where R(9) represents the 2x2 rotation matrix that rotates the neutrino state through 
an angle 9. Here 9 12 and 9 m are the mixing angles between v e and z/ 2 in vacuum 
and the mantle respectively. Clearly, #12 equals the solar neutrino mixing angle. The 
matrix $(0 m ) = diag(l, e l ^ m ) represents the change in relative phases of v\ and P 2 while 
traversing the mantle: <p m ~ 2Arr? m Ly 1 where Am^ is the mass squared difference 
between z/i and z/2 inside the mantle in units of 1CT 5 eV 2 , and L m is the distance 
traveled through the mantle in units of 1000 km. The "inverse energy" parameter is 
defined as 

y = 12.5 MeV/E (3) 

where E is the neutrino energy. The energy dependence of all quantities will always be 
implicit henceforth. 

The survival probability Eq. (J2J) may be written as 

p D « cos 2 9 12 - sin 29 m sin(20 m - 20 12 ) sin 2 (A^L m y) . (4) 

The energy dependence of p D introduces modulations in the energy spectrum of v e , 
which may be observed in the form of local peaks and valleys in the spectrum of the 
event rate aF® plotted as a function of y. The modulations are equispaced, indicating 
the presence of a single dominating frequency. These modulations can be distinguished 
from random background fluctuations that have no fixed pattern by using the Fourier 
transform of the inverse energy spectrum 

The net v e flux at the detector may be written using Eqs. and (jlj), in the form 

F° = sin 2 9 l2 F° + cos 2 12 F° + AF°A m sm 2 (k m y/2) , (5) 

where AF° = (F® — F®) depends only on the primary neutrino spectra, whereas 
A m = — sin2^ m sin(2^ m — 26*12) depends only on the mixing parameters and is 
independent of the primary spectra. The last term in Eq. (jSJ) is the Earth oscillation 
term that contains a frequency k m = 2Am^ n L m in y, the coefficient AF°A m being a 
relatively slowly varying function of y. The first two terms in Eq. (jSJ) are also slowly 
varying functions of y, and hence contain frequencies in y that are much smaller than 
k m . The dominating frequency k m is the one that appears in the modulation of the 
inverse-energy spectrum. 

The frequency k m is completely independent of the primary neutrino spectra, and 
indeed can be determined to a good accuracy from the knowledge of the solar oscillation 
parameters, the Earth matter density, and the position of the SN in the sky. 

2.3. Passage through the mantle and the core 

We now study analytically the spectral modulations arising when the neutrinos travel 
through both the mantle and the core, and the effect of the sharp density jump at their 
boundary. We denote the mixing angles, phases and mass squared differences in the core 
by replacing the superscript /subscript "m" in the last section by "c". The neutrinos 
cross two sections of the mantle with equal length L m /2 each. We denote the total 
distance traveled through the core by L c . 
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The anti-neutrino survival probability is given by 
p D = | [R(-9 m )<S>(<p m /2)R(9 m - # C )$(0 C ) 

xR(6 c -6 m )${<t> m /2)R{6 m -6 12 1 

This may be written in the form 

p D ^A + j2A l sm 2 ( ( j )l /2), 

i=l 

where A Q = cos 2 0x2 ■ The explicit expressions for the other A\ and <pi are given in 
Table [21 As in the mantle-only case, (fi m = 2Am 2 n L m y and (fi c = 2Am 2 L c y. Note that 
in the absence of travel through the core, C = and Eq. (J7J) reduces to Eq. (JU). 



(6) 



(7) 
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4>i 


M 




1 


4>m/2 


-\ sin(20 12 - A8 m ) sin(40 c - A8 m ) 




2 


(<j) m /2 + (f) c ) 


cos 2 (9 c - 6 m ) sin(26»i2 - 40 m ) sin(20 c - 26 m ) 


OH 


3 


{4>m + 4>c) 


sin(26» 12 - 26 m ) cos 4 (<9 c - 9 m ) sin(20 m ) 


OH 


4 


4>c 


- sin 2 (2# c - 29 m ) [cos(2#i 2 - 40 ro )- 
-|sin(26»i 2 -26i m )sin(26' m )] 


C(w 2 ) 


5 




| sin(26»i 2 - 26 m ) sin 2 (2(9 c - 26 m ) sin(2(9 m ) 


0(w 3 ) 


6 


(4> m /2 - <f) c ) 


-2sin(26» 12 - 40 m ) cos((9 c - 6 m ) sin 3 ((9 c - 6 m ) 


0(w 3 ) 
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{4>m - 4>c) 


sin(26> 12 - 26 m ) sm 4 (6 c - 6 m ) sin(2(9 m ) 


0(w 5 ) 



Table 2. Explicit expressions for <^ and Ai in Eq. 
The net v e flux at the detector may be written using Eqs. Q and ((Zj) in the form 



F° = sin 2 9 l2 F° + cos 2 # 12 F° + AF° £ A* sin 2 (A; i2 //2) 



(8) 



i=l 



where fcj = ^/y are the dominating frequencies. 

Not all the frequencies fcj are equally important. We estimate the relative 
magnitudes of these terms in the following manner. The mixing angles in the mantle 
and the core, 9 m and 9 C respectively, are given by 



sin 2 26 „ 



sin 2 29 



12 



(9) 



m(c) sin 2 29 12 + (cos 29 12 - 2EV m(c) / Am 2 ) 2 ' 

where V m i c ) is the effective potential due to the matter in the mantle (core) for v e . For 
the densities of the mantle as well as the core, both (9 12 — 9 m ) and (6 m — 9 C ) are small 
numbers of order 0.1 as can be seen in the left panel of Fig. In the last column of 
TableEJwe symbolically denote either of these quantities by u ~ 0.1 and show the power 
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Figure 1. Left: The energy dependence of the changes in the neutrino mixing angles 
during the passage through vacuum, mantle and core. Right: The energy dependence 
of the magnitude of A^s in Eq. (JSJ. 



of uo involved in the coefficient of that particular frequency term. Terms with higher 
powers of u are suppressed. 

We observe that the low- frequency component Ao, that does not contribute to the 
Earth effects, is the largest in magnitude. Among the terms relevant for the Earth 
effects there are three dominant frequencies, corresponding to the first three terms 
in the summation in Eq. (jHJ). The fourth term is subleading and the rest are too 
suppressed to be of any significance. The right panel of Fig. confirms that the first 
three terms have similar magnitudes, in particular \Ai\ ~ l^bl, whereas the others are 
significantly suppressed for all the energies relevant for SN neutrinos. We expect the 
first three terms to give rise to three dominant peaks in the Fourier spectrum of the 
inverse energy spectrum. Note that since Am^ ~ Am^ ~ Awig to within 20% in the 
relevant parameter range, the positions of these peaks are also known once the distance 
traversed through the Earth is known, independently of the primary neutrino spectra. 
This distance can be determined with sufficient precision even if the SN is optically 
obscured using the pointing capability of neutrino detectors [TH] . 



Signatures of supernova neutrino oscillations 



7 



3. Peaks in the power spectrum of u e 

3.1. Definitions 

The neutrino signal is observed as a discrete set of events. The measured energy of 
these events is correlated to the neutrino energy via kinematics and detector properties. 
Then we define the power spectrum of TV detected events as 



G(k) = = 



N 

i=i 



ikyi 



(10) 



In the absence of Earth effect modulations, G(k) is expected to have an average value 
of one for k > 40 ^3] • Earth effects introduce peaks in this power spectrum at specific 
frequencies, the identification of which correspond to the identification of the Earth 
effects. 

Before defining an algorithm to analyze neutrino signals from a single SN, we 
perform a check of the analytical features derived in Sec. |2 with a realistic Monte 
Carlo simulation. For the time-integrated neutrino fluxes we assume distributions of 
the form [TTi] 



F u = — 



$o (1 + a 



,1+a 



E 



£/Q. 



E r(l + a) 

where F° denotes the flux of a neutrino species emitted by the SN scaled appropriately 
to the distance traveled from the SN to Earth. Here Eq is the average energy and a a 
parameter that relates to the width of the spectrum and typically takes on values 2.5-5, 
depending on the flavor and the phase of neutrino emission. The values of the total flux 
$o and the spectral parameters a and Eq are generally different for u e , v e and u x , where 
v x stands for any of i/ MjT or v^ r . 



Model 


(Eo(Ve)) 








G 


12 


15 


18 0.8 


0.8 


L 


12 


15 


24 2.0 


1.6 



Table 3. The parameters in the neutrino spectra models motivated from SN 
simulations of the Garching (G) group and the Livermore (L) group. We assume 
a = 3 for all neutrino species. 



In order to study the model dependence, we consider two models that give very 
different predictions for the neutrino spectra. The first is motivated by the recent 
Garching calculation [20] that includes all relevant neutrino interaction rates, including 
nucleon bremsstrahlung, neutrino pair processes, weak magnetism, nucleon recoils and 
nuclear correlation effects. The second is the result from the Livermore simulation [21] 
that represents traditional predictions for flavor-dependent SN neutrino spectra that 
have been used in many previous analyses. The parameters of these models are shown 
in Table El To study the background, we use the mixing parameters of scenario B in 
Tabled in which the Earth effects are absent. 
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For the Earth density profile we use a 5-layer approximation of the PREM model 
as parametrized in Ref. ^3]. We start with a SN at 10 kpc and simulate the neutrino 
signal at a 32 kiloton scintillation detector and a megaton water Cherenkov detector. In 
both detectors, the neutrino signal is dominated by the inverse beta reaction i> e p — > ne + , 
while all other reactions have a negligible influence on the analysis below. The detector 
response is taken care of in the manner described in Refs. [TJE]. The major difference 
between the scintillation and the water Cherenkov detector is that the energy resolution 
of the scintillation detector is roughly six times better than that of the Cherenkov 
detector. This compensates the size advantage of a megaton water Cherenkov detector. 

3.2. Large scintillation detector 

To start with, we consider the power spectrum resulting from averaging 1000 SN 
simulations. This eliminates the fluctuations in the background, and illustrates the 
characteristics of the peaks in a clear manner. The power spectrum at a 32 kt 
scintillation detector for different distances traveled through the Earth is shown in Fig. El 
The top panels use the Garching model whereas the bottom ones use the Livermore 
model. The left panels show three typical cases when the neutrinos traverse only the 
mantle, whereas the right panels represent the passage of neutrinos through the mantle 
as well as the core. Only inverse beta decay events have been taken into account. 
We have checked that the inclusion of the other reactions like the elastic scattering off 
electrons or the charged current reactions on carbon do not change the results. The 
neutral current reactions on carbon, although providing a large number of events, have 
been neglected, as the monoenergetic photon produced by the decay of the excited 
carbon could be tagged in such a detector [T7] . 

The following observations may be made: 

(i) The average background approaches 1 for k > 40, as expected jTHj. The region 
k < 40 is dominated by the "0-peak," which is a manifestation of the low frequency 
terms in Eqs. (J3J) and (JHJ). Note that the 0-peak in the background case is wider 
than that in the signal case. This is because the background case corresponds to 
the scenario B, which is also the one wherein there is a complete interchange of the 
v e and v x spectra. The energy of the detected v e spectrum is thus higher, which 
results in a broader 0-peak. 

(ii) When the neutrinos traverse only the mantle, only one peak appears at the expected 
value of k m that is proportional to the distance L m traveled through the mantle. 
For L m < 3000 km, the position of the peak lies at such low frequencies that it can 
hardly be distinguished from the 0-peak. This illustrates that neutrinos must travel 
a minimum distance through the Earth before the Earth effects become observable. 

(iii) When the neutrinos travel also through the core, we observe three dominant peaks 
in each case, corresponding to k\, k% and k$ in Eq. (JHJ). We observe that, as the 
total distance traveled through the Earth increases, 
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• the third peak moves to higher k, since k 3 oc (L m + L c ). 

• the second peak, whose position is proportional to (L m /2+L c ), also gets shifted 
towards higher k as the increase in L c is larger than the decrease in L m /2. 

• the first peak, on the other hand, moves to lower k values, since as the 
trajectory of neutrinos approaches the center of the Earth, the distance traveled 
through the mantle decreases and so does the frequency of the lowest peak, 
k\ oc L m /2. This makes the detection of this peak harder at larger distances 
traveled through the Earth. 

(iv) The model independence of the peak positions may be confirmed by comparing the 
top and bottom panels of Fig. |21 The peaks obtained with the Livermore model 
are strong result of the larger difference between the v e and v x spectra in 

that model, which increases the value of AF° in Eqs. (JSJ) and (jHJ). However, the 
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Figure 2. Averaged power spectra in the case of a large scintillator detector for 
different SN models, Garching (G) and Livermore (L), and distances travelled through 
the Earth. 
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positions of the peaks are the same as those obtained with the Garching model. 
3.3. Megaton water Cherenkov detector 

The energy resolution of a water Cherenkov detector is about a factor of six worse than 
that of a scintillation detector. This means that the energy spectrum is more "smeared 
out" and higher frequencies in the spectrum are more suppressed. This makes the peak 
identification more difficult, and even a detector of the size of Super-Kamiokande turns 
out not to be sufficient [T3] . We show the power spectrum expected at a megaton water 
Cherenkov detector in Fig. El for the two SN models considered here and for different 
locations of the SN. Again we have only assumed the events from the inverse beta decay, 
since the contribution from other reactions is significantly smaller. 

The following observations may be made from the power spectra: 
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Figure 3. Same as the Fig. El but for the case of HK. 
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(i) The O-peak is broader, which makes the averaged background approach one for 
larger values of k (k c± 50) than in the case of the scintillation detector. This is 
because the energy smearing decreases the strengths of high frequency components 
and increases the strengths of low frequency components. 

(ii) In the case of neutrino propagation only through the mantle, the peak shifts to 
higher k as L m increases, as expected. However, the suppression of high frequencies 
tends to shift the peak locations to slightly lower k values than at the scintillation 
detector. In addition, as the peak position moves to larger k m , the strength of the 
peak decreases and the peak becomes harder to detect. 

(iii) When the neutrinos traverse both the mantle and the core, it is observed that 

• the third peak among the expected three dominant ones, corresponding to 
k 3 oc (L m + L c ), is highly suppressed due to large k 3 and is undetectable. 

• the other two peaks corresponding to k\ oc L m /2 and k 2 oc (L m /2 + L c ) have 
lower k values, and are not as suppressed as the k 3 peak. Moreover, the k\ 
peak is stronger than the k 2 one. 

• the k\ peak moves to lower frequencies as the total distance traveled through 
the Earth increases. Beyond a certain distance, it merges with the background 
0-peak and becomes undetectable. 

(iv) The peak positions with the Garching and Livermore spectra are at nearly the same 
positions, though the extra feature of high k suppression is observed to shift the 
peaks with the Garching model to slightly lower values of k as compared to those 
with the Livermore model. The peaks with the Livermore model naturally have 
more strength than the ones with the Garching model. 

4. Distinguishing the peaks from the background 

4-1. An algorithm for peak identification 

Though the analytic approximations seem to work well with the averaged power 
spectrum, the understanding of the statistical fluctuations within the signal from a single 
SN is crucial for the identification of the peaks. As observed in ^Hj and confirmed in 
Sec. El the average of the background power spectrum is indeed one for all values of k 
after the dominant low frequency peak in the power spectrum dies out. As long as we 
are free of the influence of this low frequency peak, the area under the power spectrum 
between two fixed frequencies k min and /c max is on an average (A; max — k m i n ). In the 
absence of Earth effects, this area will have a distribution centered around this mean. 
The Earth effect peaks tend to increase this area. If the area in a specific interval is 
found to be more than what mere background fluctuations can allow, the peak can be 
identified with confidence. 

Since in the real world the presence of fluctuations in the signal will spoil any naive 
theoretical peak, we need to introduce a prescription to carry out the analysis. 
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1. Once we know the total distance traveled by the neutrinos through the Earth, 
we can calculate the position where the peak should lie. Instead of looking for 
the maximum in the height of the power spectrum, we consider a more robust 
observable, namely the area around the position of the peak, as illustrated in Fig.|U 
When only one peak centered at k m is expected, we consider the interval k m ± Ak 
with Ak = 30, roughly the expected width of the peak. In order to avoid the 0-peak 
contamination we set a lower limit at k = 40 (k = 50) for the scintillator (water 
Cherenkov) detector. When the neutrinos also cross the Earth core, multiple peaks 
are present and we measure the area from k = 40 until k = 160 in such cases. 

2. The next step is to analyze the statistical significance of the result obtained. For 
this purpose we compare the value of the measured area with the distribution of 
the area in the case of no Earth matter effects. Since the different frequencies are 
not uncorrelated, the background distribution is not simply Gaussian centered at 
A av = (&max — ^min) and with width V^max — ^min- Therefore, we perform a Monte 
Carlo analysis of the background case and calculate the exact distribution with 
which one can compare the actual area measured. We illustrate this in the right 
panel of Fig. 01 where the black curve shows the area distribution of the background. 
The confidence level of peak identification may then be defined as the fraction of the 
area of the background distribution that is less than the actual area measured. We 
denote the area corresponding to a% C.L. by A a . Figure |U also shows A 95 « 100, 
the area corresponding to a peak identification with 95% confidence. 




Figure 4. Left: Realistic spectrum from a single simulation. Right: Area distribution 
of the background (black) and the signal (red) obtained for a 32 kton scintillator 
detector and Garching model for r\ = 60. 
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4-2. Quantifying the efficiency of the algorithm 

Since the peak identification algorithm is statistical in nature, it is worthwhile to have 
an idea of the probability with which a peak can be identified with a given confidence. 
This probability clearly depends on the distance traveled by the neutrinos through the 
Earth, which in turn is determined by the location of the SN in the sky. We parameterize 
the SN location by the nadir angle rj of the SN direction at the detector. 

We simulate the area distribution for the "signal" using the neutrino mixing 
scenarios that allow Earth effects and compare it with the background distribution. 
The probability p a of peak identification at a% CL. is the fraction of the area of the 
signal distribution above A a . In Fig.QJ it has been indicated with the red-hashed region. 
In Fig. El we show p 95 and p 99 as a function of rj in the case of a scintillation detector, for 
the two SN models considered. An increase in the distance traveled through the Earth 
corresponds to a decrease in rj. The passage through the core corresponds to r] < 33°. 
One can see that the presence of the core enhances the chances of detecting the Earth 
matter effects. The probability is higher at values of rj close to the boundary between 
the mantle and the core because the three peaks are clearly visible. The oscillation 
pattern arising at this region stems from the interference of the first two peaks in the 
spectrum, whose positions at this point differ only by C oc L c . As the distance traveled 
by the neutrinos through the core increases, p a decreases, the reason being the approach 
of the first peak to the lower limit frequency, k ~ k min , and its eventual disappearance. 




Figure 5. Left: Comparison of pg$ and pgg for the Garching (G) and Livermore (L) 
SN models in a 32 kton scintillator detector. Right: Comparison of pg$ in this large 
scintillator detector (SC) and in the case of a megaton water Cherenkov (HK), for the 
Garching model. 
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As expected, the chances of peak identification are also higher when the primary 
spectra of v e and v x differ more. As shown in Fig. El the Livermore model predicts much 
larger chances of a successful peak identification. 

In the right panel of Fig. we assume the Garching model and compare the results 
obtained with a 32 kton scintillator detector and a megaton water Cherenkov detector. 
In the latter neutrinos travel more and more distance in the mantle the peak 

moves to higher k values, and due to the high k suppression as described in Sec. 13.31 the 
efficiency of peak identification decreases. However, when the neutrinos start traversing 
the core, additional low k peaks are generated and the efficiency increases again. 

One of the features of this algorithm is its robustness. However in some cases it 
turns out to be very conservative. For instance, if we have a look at the megaton water 
Cherenkov, the peaks at the lowest frequencies are almost eaten-up by the 0-peak of the 
background when the neutrinos cross the Earth core, cf. Fig. 01 Setting A; min = 50 as 
the lower limit of the area integration results therefore in a loss of considerable amount 
of information under these conditions. In this particular case it is possible to optimize 
the efficiency of the method by choosing a floating lower cut: instead of considering a 
fixed value, k m i n = 50, as the lower limit for the area integration, one defines k m [ n as the 
frequency at which the spectrum has the first minimum after neglecting the effect of 
spurious fluctuations. With this modified prescription one can again compare the area 
distribution for the background and that of the signal, and calculate a new p a . We have 
checked this method for both the scintillator and the water Cherenkov detector. For 
the former the improvement is not relevant. The reason is that the second and third 
peaks contribute significantly to the signal. So, even when the first peak disappears 
behind the 0-peak of the background for fixed lower cut-off, the loss of information 
is not important. However in the case of a water Cherenkov detector the efficiency 
is significantly enhanced for paths traversing the Earth core as can be seen in Fig. |U1 
Due to the suppression of the peaks located at higher frequencies only the first peak 
contributes significantly to the signal. However for trajectories involving small rj this 
peak is centered at very low frequencies, almost completely hidden by the 0-peak of 
the background. Under this situation if one allows the lower bound to shift to smaller 
frequencies the whole peak contributes to the signal, and therefore the probability to 
see the Earth matter effects increases. On the other hand when the neutrinos only cross 
the mantle the location of the unique peak is mostly at k > k min = 50. Therefore, 
this modified prescription does not help much to improve the efficiency to observe the 
modulation of the neutrino spectra due to the Earth matter effects. 

5. Summary and conclusions 

When neutrinos coming from a core-collapse supernova pass through the Earth before 
arriving at the detector, the spectra may get modified due to the Earth matter effects. 
The presence or absence of these effects can distinguish between different neutrino 
mixing scenarios. We have seen that these Earth matter effects on supernova neutrinos 



Signatures of supernova neutrino oscillations 



15 




5 10 15 20 25 30 



7] 

Figure 6. Comparsion of pg$ for two different methods of integration, with fixed 
and floating lower limit, in the case of a megaton water Cherenkov detector and the 
Garching model. 

can be identified at a single detector through peaks in the Fourier transform of their 
"inverse energy" spectrum. The position of these peaks are independent of the initial 
neutrino fluxes and spectral shapes. 

We have performed an analytical study of the positions and the strengths of the 
frequencies that characterize the inverse-energy spectrum of the neutrinos for different 
neutrino trajectories through the Earth. In the case that the SN neutrinos only traverse 
the mantle a single peak shows up in the power spectrum. In contrast we have observed 
that if both the mantle and the core are crossed before the neutrinos reach the detector as 
many as seven distinct frequencies are present in the inverse energy spectrum. However 
only three peaks are dominant in the power spectrum. This increase in the number of 
expected peaks leads to an easier identification of the Earth matter effects. 

In order to illustrate the qualitative features of the present analysis we have 
considered the power spectrum resulting from averaging 1000 SN simulations for 
different SN models and different detector capabilities. In particular we have assumed a 
32 kton scintillator detector and a megaton water Cherenkov detector. We have shown 
how the energy resolution turns out to be crucial in detecting the modulation introduced 
in the neutrino spectra by the Earth matter effects. First, the better resolution of the 
scintillator detector compensates for the larger water Cherenkov detector size. On the 
other hand, the worse energy resolution in water Cherenkov detectors does not only 
imply the need of a larger volume but also suppresses significantly the peaks at higher 
frequencies, in contrast to the case of scintillator detectors. 

We have considered two different SN models as an illustration of the current 
uncertainties in the initial fluxes. We have observed that the strength of the peaks is 
larger in those SN models with bigger differences between v e and spectra. However, 
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we have found that the position of the peaks is model independent. Therefore their 
identification serves as a clear signature of the Earth matter effects on SN neutrinos, 
which in turn can help to discard the neutrino mass scheme with inverted mass hierarchy 
and sin 2 G 13 > KT 3 . 

We have introduced a simple algorithm to identify the peaks in the presence of 
background fluctuations. This method is based on the integration of the area around 
the expected position of the peak. By comparing the area distribution without and 
with the spectral modulations induced by the Earth matter effects we have analyzed 
the statistical significance of the result. As expected the presence of the core as well 
as a larger difference in the initial spectra enhance the probability of identifying the 
Earth effects. We have also presented a variation of the algorithm which improves its 
efficiency significantly in the case of a water Cherenkov detector for neutrino trajectories 
passing through the core. Therefore we believe that more efficient algorithms could be 
developed. The results we presented should be considered therefore as conservative 
lower limits. 
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